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We search for invisible decays of the T(1S) meson using a sample of 91.4 x 10® T(35) mesons col- 
lected at the BAa4R/PEP-II B Factory. We select events containing the decay r(3S') 7r+7r"T(15') 
and search for evidence of an undetectable T{1S) decay recoiling against the dipion system. We set 
an upper limit on the branching fraction B{T{1S) — > invisible) < 3.0 x 10~* at the 90% confidence 
level. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The nature of dark matter is one of the most challeng- ing issues facing physics. Observation of standard model 
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(SM) particles coupling to undetectable (invisible) final 
states might provide information on candidate dark mat- 
ter constituents. In the SM, invisible decays of the T{1S) 
meson proceed by bb annihilation into a vV pair, with a 
branching fraction B{T{1S) invi sible) « 1 X 10-^ [H, 
well below the current experimental sensitivity. However, 
low-mass dark matter candidates could couple weakly to 
SM particles to enhance the invisible branching fraction 
to the level of IQ-* to lO^^ [2]. 

Searches for this decay of the T{1S) can be carried out 
at e+e~ colliders operating at the T{2S) or T{3S) reso- 
nance. The presence of the T{1S) is tagged by detecting 
the particles emitted in decays of the resonance to T{1S). 
Previous searches by the CLEO Q and Belle Q collab- 
orations yielded upper Hmits of B{T{1S) invisible) < 
3.9 X 10"^ and < 2.5 x 10"^ at the 90% confidence level 
(C.L.), respectively. In this paper we present a search for 
this final state using almost an order-of-magnitude more 
T{1S) mesons. 

The data used in this analysis were collected with 
the BaBAR detector at the PEP-II asymmetric-energy 
e~^e~ collider running at an e"*"e~ center-of-mass 
(CM) energy corresponding to the mass of the T{3S) 
(10.3552 GeV/c2 d). The presence of a T(IS') meson is 
tagged by reconstructing the tt+tt" pair (dipion) in the 
transition 7(35) TT+Tr-r{lS). The BaBAR detector is 
described in detail elsewhere [6|. These data were taken 
using an upgraded muon system, instrumented with both 
resistive plate chambers |6| and limited streamer tubes 
between steel absorbers [7|. For these data the trigger 
was modified to substantially increase the dipion trigger 
efficiency for the signal process. The data sample con- 
taining these improvements represents 96.5 x 10^ T{3S) 
mesons. 

We model both generic T{3S) decays and the signal 
process using a Monte Carlo (MC) simulation based on 
Geant4 [a]. The r{3S) -> 7r+7r-r(lS') events are gen- 
erated according to the matrix elements measured by the 
CLEO collaboration In signal events, the mass recoil- 
ing against the dipion (Mrec) peaks at the T{1S) mass 
(9.4603 GeV/c^ [H]). The same is true for background 
events in which a real T{3S) 7r+7r^T(lS') transition 
occurs but the T{1S) final-state particles are undetected 
("peaking background"). However, the dominant back- 
ground containing a pair of low-momentum pions does 
not exhibit this structure ("non-peaking background"). 
The analysis strategy is as follows: first apply selec- 
tion criteria to suppress background, primarily the non- 
peaking component; then fit the resulting Mj-oc spectrum 
to measure the peaking component (signal plus peaking 
background) . 

We define three subsamples for both data and MC 
events. The first of these, the "invisible" subsample, is 
designed to contain signal events. Events in this sub- 
sample have only two charged pions. The other two sub- 
samples are used to check and correct MC predictions 



for the processes which contribute to the peaking back- 
ground in the invisible subsample. The "four-track" sub- 
sample consists of events with two pions plus two tracks 
with high momenta in the CM frame, consistent with 
two-body decay of the T{1S) where both final-state par- 
ticles are detected. The "three-track" subsample con- 
sists of events containing two pions and only one high- 
momentum track, consistent with two-body T{IS) decay 
where only one of the final-state particles is detected. 

We select events in the invisible subsample by requir- 
ing that there are exactly two tracks originating from 
the interaction point ("IP tracks"), with opposite elec- 
tric charge. An IP track is required to have a point of 
closest approach to the interaction point within 1.5 cm 
in the plane transverse to the beams and within 2.5 cm 
along the z-axis. We further require these tracks to 
each have CM-frame momentum p* < 0.8GeV/c, con- 
sistent with pions from the dipion transition. The dip- 
ion system is required to have an invariant mass satisfy- 
ing Mtttt G [0.25, 0.95] GeV/c^, compatible with kinematic 
boundaries (M^^ e [2M^, {Mr{3S) - Mr{is))]) after al- 
lowing for resolution effects. The dipion recoil mass is 

M^,, = s + M^^-2V^E:^, (1) 

where E*^ is the CM energy of the dipion system 
and = 10.3552 GeV/c^. We require that M,ec e 
[9.41,9.52] GeV/c^. The efficiency of this selection for 
signal events is about 64%, due to the requirement of 
reconstructing the two pions. All selection criteria were 
finalized without looking at data in a narrower M^ec "sig- 
nal region" that, according to simulation, contains more 
than 99% of the signal (see discussion of signal shape 
below for the precise signal-region definition). 

We select three-track and four-track events using the 
same dipion selection as in the invisible subsample. We 
search for high-momentum tracks from the T(1S) decay 
(i.e., from T{1S) e+e" or T(IS') We re- 

quire that there be one or two additional IP tracks, each 
with p* > 2.0GeV/c. If either of these tracks passes 
electron identification criteria, both are treated as elec- 
trons; otherwise, both are treated as muons. In the for- 
mer case we account for possible radiative energy loss 
due to bremsstrahlung by pairing an electron with a pho- 
ton emitted close in angle and increasing the electron's 
energy and momentum by the energy of this photon. 
When two high-momentum tracks are present, we require 
that they have opposite charge and a two-track invariant 
mass €E [9.00, 9.80] GeV/c^. We remove photon conver- 
sions from these events by rejecting the event if either 
pion satisfies electron-identification criteria. This intro- 
duces a negligible efficiency loss: the probability of a pion 
to be misidentified as an electron is « 0.1%. Finally, we 
require that the mass difference between the 7r"'"7r~£"'"^~ 
and t+t- systems e [0.89,0.92] GeV/c^. 

At this stage, the background level in the invisible sub- 
sample is several orders of magnitude larger than any 
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hypothetical signal. We reject most of this remaining 
background with a multivariate analysis (MVA), imple- 
mented as a random forest of decision trees [10|. The 
random forest algorithm is trained on signal MC events 
and 5.3% of data outside the signal region in M-^cc- The 
contribution of peaking components to these data is neg- 
ligible. The data and signal MC events used to train the 
MVA are excluded from the rest of the analysis, leaving 
91.4 X 10^ T{iS) events for use in the final result. 

We use the following variables, which have been deter- 
mined to be only weakly correlated with Mpec, as inputs 
to the MVA: the probability that the pions originate from 
a common vertex; the laboratory polar angle and trans- 
verse momentum of the dipion system; the total number 
of charged tracks, IP tracks or otherwise, reconstructed 
in the event; booleans that indicate whether either pion 
has passed electron, kaon, or muon identification crite- 
ria; the cosine of the angle (in the CM frame) between 
the highest-energy photon (71) and the normal to the 
decay plane of the dipion system; the energy in the lab- 
oratory frame of the 71; the total neutral energy in the 
CM frame; and the multiplicity of candidates, defined 
using the shape and magnitude of the shower resulting 
from interactions in the calorimeter. 

The selection on the MVA output is optimized by 
choosing the threshold that achieves the minimum 
statistical uncertainty (dominated by background) on 
B{T{1S) invisible) and, in the null signal hypothe- 
sis, the lowest upper limit at the 90% C.L. Both were 
achieved by requiring an MVA output > 0.8 (where the 
full range is to 1). The efhciency of this criterion for 
signal-MC events is 37.0%, as compared to 0.8% for data 
events outside the signal region. The total efficiency of 
all trigger and event selection requirements is determined 
from signal-MC simulation to be 16.4%. 

Figure 1 shows the Mrec distribution for the selected 
events. We extract the peaking yield by an extended 
unbinned maximum likelihood fit, with the non-peaking 
background described by a first-order polynomial. The 
signal and the peaking background should have the same 
shape in Mroc- We describe this shape using a mod- 
ified Gaussian function with a common peak position 
(/io), independent left and right widths {aL,R), and non- 
Gaussian tails (governed by parameters aL,Fi.)- The func- 
tional form on either side of the peak is 

exp[-(Mrcc - iJLo?/{2al j^ + aLM{Mrcc - Mo)^)] • (2) 

We determine the parameters of this probability density 
function (PDF) by fitting Afrcc in the four-track data sub- 
sample. The signal region, excluded when training the 
MVA, is defined as the region in Fig. [T] which is < ^(Jl,r 
from the peak position, Mrec S [9.4487,9.4765] GeV/c^. 

The fit to the invisible subsample then determines all 
of the parameters of the non-peaking background PDF, 



the yield of the non-peaking background, and the yield 
of the peaking component. The result for the peaking 
yield is 2326 ± 105 events. 

Using a second-order polynomial for the non-peaking 
background results in no change in the extracted peaking 
yield. The systematic uncertainty on that yield associ- 
ated with the fixed parameters in the signal PDF is es- 
timated by varying those parameters in the fit. We find 
an uncertainty of 18 events. 
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FIG. 1: The maximum likelihood fit to the dipion recoil mass 
for data in the invisible subsample. The components of the 
fit are the non-peaking background (dash-dotted line) and 
the peaking component (solid filled). The total fit function is 
also shown (solid line) . Comparing the fit to the binned data, 
X^/(degree of freedom) = 0.973. 

We next estimate the contribution of background to 
the peak. The MC simulation predicts 1019 T(15) 
e+e- events, 1007 Y{IS) events, 92 T(15) ^ 

T+r^ events, and 2.9 ± 1.3 T{\S) hadrons events. 
These predictions depend upon branching fractions 
which have significant uncertainties Q and on the ac- 
curacy of the modeling of event reconstruction and se- 
lection. We use four-track and three-track data and MC 
subsamples to test and correct the MC prediction of 2122 
total events. 

We first use the four-track subsamples to calibrate the 
product of the branching fractions for T{1S) — > £~^i~ and 
T{3S) 7r+7r~T(lS') and the dipion reconstruction effi- 
ciency. We compare the event yields between four-track 
data and MC subsamples when the positively-charged 
lepton is emitted in the central section of the tracking 
system, |cos(6';+)| < 0.3 (laboratory-frame angle). The 
simulation underestimates the number of events in data 
by a factor of (1.088 ± 0.012). This is plausible in light 
of the branching fraction uncertainties and track recon- 
struction uncertainties (« 0.5% per track). Since the 
effect of the high-momentum track reconstruction has a 
negligible contribution here, this data/MC correction fac- 
tor is applied to all of our MC-simulation subsamples. 
For the four-track subsample, Fig. UK a) shows that the 
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distribution of the high-momentum tracks in the detec- 
tor is very well described by the MC simulation at all 
polar angles. 

We next compare the data and MC efficiencies for re- 
constructing the single lepton in the three-track subsam- 
ple. Any discrepancy would imply a complementary mis- 
take in the invisible peaking background. Given the CM- 
frame polar angle coverage of the detector, for three-track 
events the high-momentum lepton in the forward direc- 
tion often escapes detection and thus the detected lepton 
is in the backward direction. We compare the MC and 
data laboratory-frame polar angle distributions for these 
events in Fig. [2^b). The three-track subsample, in con- 
trast to the four-track subsample, has a significant non- 
peaking background in recoil mass. Hence three-track 
peaking yields vs. polar angle are determined by using the 
Mj-cc fit described above and applying an event-weighting 
technique jlH. The MC simulation describes the distri- 
bution well everywhere except at cos(6'£) < —0.84, where 
the simulation overestimates the reconstruction rate. 

For leptons in this far-backward region, we use the ra- 
tio of data to simulation vs. lepton cos{9) from Fig. ^h) 
as the basis of an accept-reject method applied to the 
high-momentum track. When this method removes the 
track, it in effect reassigns a three-track event to the in- 
visible category. We also weight the reassigned events by 
the ratio of simulated trigger efficiencies for the three- 
track and invisible subsamples and assign 100% uncer- 
tainty to this difference in trigger efficiency. Applying 
this additive correction after the scaling correction (from 
the four-track subsample), the total peaking background 
estimate increases from 2122 events to (2451 ±38) events. 
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FIG. 2; The distribution of cos(S), where 9 is the laboratory- 
frame polar angle of (a) the positively charged high- 
momentum track in the four-track subsample and (b) the sin- 
gle high-momentum track in the three-track subsample. The 
normalization correction from the four-track subsample has 
been applied to the MC yields in both cases. 

We test the prediction of the contribution of non- 
leptonic T{IS) decays to the peaking background using 
an additional control sample. Events in this sample con- 
tain only two tracks (the pions) and pass all other criteria 
for the invisible subsample, except that the MVA require- 
ment is replaced by a requirement that the 71 has energy 
> 0.250 GeV. This selects a set of events which is almost 
orthogonal to the signal selection, since the MVA > 0.8 



requirement results in a steep falloff in efficiency vs. 71 
energy near 0.250 GeV. We compare this energy distribu- 
tion in data (using the weighting technique 11|) to that 
from simulation. As the 71 energy approaches 0.250 GeV 
from above, we find the MC simulation underestimates 
the data by no more than a factor of four. Since the ex- 
pected contribution of these events to the peaking back- 
ground is 0.14% of the total, we assign 0.6% (15 events) 
as an additional systematic uncertainty on the peaking 
background, for a total of ±41 events. 

A number of multiplicative systematic corrections and 
uncertainties to the peaking background also enter, in a 
fully-correlated manner, when the extracted signal yield 
is converted to the T{1S) invisible branching fraction. 
The ffist such contribution is the 1.088±0.012 correction 
factor derived from the four-track subsample. But this 
does not account for trigger and MVA effects which might 
differ for the invisible and four-track subsamples. Since 
events used to train the MVA have already passed the 
trigger requirements, we first study the effect of trigger 
selection on data. The EaBAR trigger consists of a hard- 
ware and a software stage. The latter is tested by using 
a heavily-prescaled sample of events which bypassed it. 
We apply the software- level trigger to these events and 
find that the ratio of efficiencies in data and MC simula- 
tion is 0.997 ± 0.009. This ratio is taken as a correction 
to the signal efficiency and the peaking background. To 
assess how well the impact of the hardware trigger on the 
two pions is simulated, four-track events are used, since 
their trigger decision is based largely on the two high- 
momentum lepton tracks. We apply to the pions a set of 
selection criteria which approximate those applied by the 
hardware trigger. The data and MC efficiencies for these 
requirements differ by 2.2%. Since this test is done on 
samples for which the hardware trigger is only approxi- 
mated, we take this difference as a systematic uncertainty 
rather than apply a correction for it. 

After applying the approximate hardware trigger cri- 
teria to the four-track subsamples for both T{3S) MC 
simulation and data, we apply the nominal MVA selec- 
tion to both. The relative difference in efficiency between 
these MC and data subsamples is 4.0%. Since the hard- 
ware trigger is again only approximated for this test, we 
apply no correction for the difference, but assign it as a 
systematic uncertainty on the MVA selection. 

Adding the multiplicative uncertainties in quadrature, 
the total correlated systematic uncertainty is 4.8%. The 
final corrected prediction for the peaking background is 
(2444 ± 123) events, including the uncorrelated uncer- 
tainty of 41 events. From this we determine the sig- 
nal yield to be (—118 ± 105 ± 124) events, where the 
errors are statistical and systematic, respectively. To 
obtain B{T{1S) invisible) we divide this by the sig- 
nal efficiency, the number of T{3S) mesons, the branch- 
ing fraction for the dipion transition (4.48% [H']) and the 
correction factors (1.088 x 0.997). The factor derived 
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from the four-track subsample includes a possible ad- 
justment of B{T{1S) e+e") + B{r{lS) Ai+M"), 
not relevant for signal. We take this adjustment to be 
1.000 ±0.025 and remove it by assigning an additional 
systematic uncertainty of 2.5%. Taking correlations into 
account, we determine that B{T{1S) invisible) = 
(-1.6± 1.4(stat) ± 1.6(syst)) x lO"'^. 

Lacking evidence for this decay, we use a Bayesian 
technique to set an upper limit on the branching frac- 
tion. We convolute the statistical likelihood, a function 
of B(T(1S) invisible), with Gaussian functions repre- 
senting the systematic error. We assume a prior proba- 
bility that is flat in branching fraction and integrate the 
likelihood from to a value such that 90% of the to- 
tal integral above is enclosed. The resulting limit is 
B(r{lS) invisible) < 3.0 x lO^'' at the 90% C.L. 

In conclusion, we search for invisible decays of the 
Y{1S) meson. We do so by looking for evidence of the 
decay of the T{1S) into undetectable final states recoil- 
ing against the dipion system in T(35) — *■ 7r+7r~T(lS'), 
using a sample of 91.4 x 10^ T{3S) mesons. We find no 
evidence for T{IS) — > invisible and set an upper limit on 
its branching fraction at 3.0 x 10"'* at the 90% C.L. This 
limit is almost an order of magnitude closer to the SM 
prediction than the best previous limit. 
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